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ABSTRACT 

We discuss the observable consequences for the detection of galaxies in the X-ray bandpass resulting 
from a peak in the cosmic star-formation rate at a redshift > 1. Following White & Chosh, we assume 
a large evolution in the X-ray/B luminosity ratio at z ~ 0.5 — 1.5 resulting from the X-ray binaries that 
have evolved from stars formed at z > 1 — 2. Using the HDF-N redshift survey data and the locally 
observed X-ray/B luminosity ratio as a guide, we estimate a median X-ray flux (2-10 keV) on the order 
of 8 X 10^^^ erg cm^^ s^^ for galaxies in the HDF-N, which is consistent with a signal derived from 
a stacking analysis of the HDF-N Chandra data by Brandt et al. (2001). We also predict the number 
counts in deep X-ray surveys expected from normal galaxies at high redshift. 

Subject headings: evolution-galaxies:evolution-X-rays: galaxies-X-rays 



1. INTRODUCTION 

The X-ray emission of nearby normal spiral galaxies has 
been observed to be dominated by X-ray binaries (XRB; 
see Fabbiano 1989 for a review). There is a strong correla- 
tion between the optical and X-ray luminosities of galaxies 
which is also consistent with an XRB population dominat- 
ing the X-ray emission (i.e., assuming that XRB constitute 
an constant fraction of stars in a galaxy). Deep opti- 
cal, IR and sub-mm surveys suggest that the peak of the 
cosmic star-formation lies at redshifts on the order of 1- 
2 or more (Blain et al. 1999). White & Ghosh (1998, 
hereafter WG98) predicted that X-ray binaries that form 
from this early star-formation activity should result in the 
enhancement of X-ray emission in galaxies, most notably 
around a redshift of 0.5-1.0 as the massive stars super- 
nova forming low-mass X-ray binaries (LMXRBs) contain- 
ing neutron stars, or blackholes The LMXRBs contain- 
ing neutron stars in turn evolve to form millisecond radio 
pulsars (MRP), resulting in the observed distribution of 
XRB and MRP in the Calaxy and nearby galaxies. Fur- 
ther modeling by Ghosh & White (2001; hereafter GWOl) 
has been undertaken to refine their previous results, with 
the primary improvement being the use of more recent 
star- formation rate (SFR) history estimates (e.g., Blain 
et al. 1999). These more recent SFR history estimates 
take into account sub-mm data which is sensitive to dusty 
(and hence obscured) high-redshift galaxies that would be 
missed by optical/UV surveys. The main impact of these 
modifications is that the current estimate of the enhance- 
ment of X-ray flux from galaxies (due to LMXRB) peaks 
at z ^ 1.5 rather than z ~ 1.0. The modeling of WG98 
and GWOl also includes the predictions for enhancements 
of high-mass X-ray binaries (HMXRB). We defer evalu- 



ation of the HMXRB enhancement for future work, but 
note here that since HMXRB have a much shorter evolu- 
tionary time scale than LMXRB, the HMXRB enhance- 
ment (along with other effects of starburst activity) would 
closely track the SFR history and therefore occur at higher 
redshifts than the LMXRB enhancement discussed here. 

A test of these X-ray evolution models would be to 
evolve the local universe X-ray luminosity function and 
derive the number counts of galaxies expected at given 
flux levels (i.e., in a logN-logS diagram) for comparison 
with deep X-ray surveys. Unfortunately, the X-ray lumi- 
nosity function of nearby galaxies is not known directly 
because the typical X-ray flux of galaxies is of order of 

logFx = —13 14 ergs cm~^ s~^ which is below the 

limiting flux of existing X-ray all-sky surveys (e.g, the 
ROSAT All-Sky Survey has a limiting flux on the order 

of logFo.5_2.okcV ~ -12 13 ergs cm"^ s~^; Voges et 

al. 1999). Here we take an alternative approach which is to 
use the known X-ray/optical luminosity correlation given 
in David, Jones & Forman (1992), determined from Ein- 
stein data, to estimate the X-ray luminosities and fluxes 
of galaxies in the HDF-N, where optical luminosities and 
redshifts have been established. Since the expected X-ray 
flux of high-redshift galaxies is low, Brandt et al. (2001) 
performed a stacking analysis of a Chandra observation of 
the HDF-N. We show below that the X-ray flux distribu- 
tion derived from our analysis is consistent with the galaxy 
flux signal detected by Chandra. A secondary goal of this 
paper is to estimate the logN-logS distribution of galaxies 
(based on the HDF-N) in order to predict X-ray number 
counts that could be detected by more sensitive surveys. 



2. METHODOLOGY 



^ Present address: The Johns Hopkins University, Dept. of Physics and Astronomy, 3400 N. Charles St., Baltimore, MD 21218; 
ptak@pha.jhu.edu 



2 



As discussed above, our main goal is to estimate the X- 
ray flux and luminosity for galaxies based on their optical 
luminosity, the X-ray/optical luminosity correlation and 
X-ray luminosity evolution models. The HDF-N "proper" 
sample presented in Cohen et al. (2000) contains rcdshifts 
and R-magnitudes for 125 galaxies spanning 4.75 arcmin^. 
The X-ray/optical correlation given in David, Jones & 
Forman (1992) is based on B-band luminosities, so ac- 
cordingly the observed R-band magnitudes must be con- 
verted to rest-frame B-band magnitudes. We used the k- 
correction plots in Frei & Gunn (1994) which include R to 
B band corrections based on optical and UV observations 
of galaxies (which we caution only extends to z=0.6), from 
which we derived the B-band luminosity for each galaxy 
(we used Hq = 70 km Mpc~-^ and qo = 0.1 through- 
out this paper). Here we are primarily interested in the 
mean properties of galaxies and accordingly did not seg- 
regate by galaxy type, which we differ to future work (sec 
also below). The k-corrcctions as a function of galaxy type 
vary by ~ 1 magnitude, which is somewhat larger than the 
scatter found for control galaxies by Frei & Gunn (1994), 
and accordingly we assume an uncertainty of 1 magnitude. 

These B luminosities were then converted to X-ray lumi- 
nosities using the Lq. 5-4.5 kev/LB relation given in David, 
Jones & Forman (1992). Assuming a power-law spec- 
trum with a photon index of 1.8 and neutral absorption 
Nh < 10^^ cm~^ (the galactic column towards the HDF- 
N is 1.62"cm"2; Stark ct al. 1992), which is consistent 
with an XRB-dominated spectrum and the observed 2-10 
keV X-ray spectra of nearby galaxies (c.f., Ptak et al. 1999 
and references therein), -F2-10 koV = -Fo.5-4.5 kcV to within 
~ 15%. Similarly, these same spectral assumptions also 
imply a F0.5_2.0koV flux that is (0.5 - 0.6)^0.5-4.5 keV 
(at redshifts exceeding 0.5, Fq. 5-2.0 keV = O.55F0.5-4.5 keV 
for all Nh < 10^-^ cm~^) and we adopt 1.0 and 0.55 as 
the conversion factor to the 2.0-10.0 kcV and 0.5-2.0 kcV 
bandpasses. We also note that the stellar evolution models 
we are concerned with here are relevant only on long look- 
back time scales, and therefore only the lower-mass, older 
stellar populations have a significant impact in our anal- 
ysis. Ignoring the younger stellar populations biases our 
results in the sense of underestimating the expected X-ray 
flux for starburst systems which tend to exhibit enhanced 
supernovae and high-mass X-ray binary populations (see 
Fabbiano 1989). Since the David, Jones & Forman (1992) 
correlation was performed using all galaxies types, the dif- 
fering amount of X-ray production as a function of galaxy 
type inherently contributes to the scatter of the correla- 
tion. Accordingly, the results of this paper should be ac- 
curate in a statistical sense but should not be applied to 
any individual galaxy. 

The resultant B and X-ray luminosity distributions de- 
rived from this procedure are shown in Figure 1. The 
errors given in this figure are derived from a Monte-Carlo 
approach in which the conversion to the B and X-ray band- 
passes was repeated (1000 iterations) with gaussian de- 
viates with a standard deviation of 1 magnitude added 
to the resulant B magnitudes and gaussian deviates with 
standard deviation of 1 added to the log X-ray luminosities 
(i.e., simulating an order-of- magnitude scatter in Lx/ Lb)- 
Errors derived in this way were comparable to the count- 
ing errors that would be expected in each histogram bin 



(i.e., vW, where N = the number of galaxies in a given 
histogram bin). The mean galaxy L2-10 keV is 3.8 x 10^® 
erg s~^, which is comparable to the mean value for pas- 
sive galaxies given in Georgantopoulos, Basilakos, & Plio- 
nis (1999) of 3.2 x 10^*^ erg s^^ (after adjusting to our 
values of Hq and qo). These values are somewhat lower 
than other nearby galaxies X-ray luminosity estimates 
(e.g., Fabbiano, Trinchieri, & McDonal 1984) however it 
should noted that here and in Georgantopoulos, Basilakos, 
& Plionis (1999) active (including narrow-line) and star- 
burst galaxies have not been (explicitly) included in the 
analyses. 

2.1. Evolution in Lx/Lb 

We proceed now assuming that the ratio Lx/Lb evolves 
as a result of an excess of X-ray binaries at earlier epochs. 
We use the evolutionary models "Peak-M" and "Gaussian" 
given in Ghosh & White (2001), plotted in Figure 2. Both 
models arc based on fits to the observed star-formation 
rate using optical/UV observations, however the "Gaus- 
sian" model also includes a component with a Gaussian 
functional form that takes into account the evolution of 
the IR luminosity function using IR and sub- mm data (c.f., 
Blain et al. 1999). These models only take into accoimt 
the relative change in X-ray luminosity as a function of 
redshift, however in this paper we are attempting to de- 
rive the X-ray luminosities of galaxies based on the rest- 
frame B-band luminosity. Accordingly, any optical evolu- 
tion must be taken into account explicitly before we could 
apply our X-ray evolution models. To this end, we took the 
mean and standard deviation of the B-band luminosities 
in redshift bins large enough to contain at least 20 galax- 
ies and fit log Lb as a function of z with a linear model 
(also shown in Figure 2). The X-ray luminosities were 
evolved by multiplying them by the ratio Ex{z) / Eb{z), 
where Ex{z) and Eb{z) are the X-ray and optical (after 
normalizing to z=0) evolution models shown in Figure 2. 
The net amount of optical evolution derived in this way 
amounts to only a factor of 2 from a redshift of to 1. 

3. RESULTS 

Applying these procedures, the mean 2-10 keV X-ray 
galaxy flux increased from 2.3 x 10~^® erg cm~^ s~^ to 

4.2 and 8.4 x 10~^^ erg cm~^ s~^ and the mean 2-10 
keV X-ray luminosity increased from 3.8 x lO''^ erg s"'^ to 

5.3 and 19. x lO^^ erg s^^ for the "Peak-M" and "Gaus- 
sian" models, respectively (with 0.5-2.0 keV values a fac- 
tor of 0.55 lower as discussed above). The distribution 
of Lx and Fx (before and after including evolution) is 
shown in Figures 1 and 3. The flux density implied is 4.3 
and 7.8 x lO^^'^ ergs s"^ cm"^ deg."^ ^j^g 0.5-2.0 keV 
and 2-10 keV bandpasses, or 2-4% of the 2-10 keV X-ray 
background. Kuntz, Snowden & Mushotzky (2001) find 
that the extra-galactic X-ray backround in the 1-2 keV 
bandpass is ^ 5.7 x 10^^^ ergs s~^ cm~^ deg.~^. With 
same spectral assumptions discussed above, Ji-2 keV ^ 
0.3i^o.5-4.5 keV) implying a 1-2 keV flux density of 
2.3 X 10"^"^ ergs s~^ cm^^ deg.^^, or ^ 4% of the observed 
1-2 keV background. We repeated our analysis after only 
including galaxies with redshifts in the range of 0.5-1.0 (55 
galaxies) which resulted in a mean 2-10 keV X-ray flux of 

7.0 X 10~^* erg cm~^ s~^ (indicating that our results are 
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Fig. 1. — The estimated B-band (left) and X-ray luminosities (right) for the HDF-N galaxies. The error bars shown are derived from 
simulations that incorporate a 1 magnitude uncertainty in the k-corrections and an order of magnitude uncertainty in the X-ray /optical flux 
ratio. The X-ray luminosity distribution is given for the unevolved case (solid line and points) as well as that derived from the "PeaJ£-M" 
(dashed line) and "Gaussian" (dotted line) evolution models. 
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Fig. 2. — (left) Evolution models "Peak-M" (solid line) and "Gaussian" (dashed line) for the X-ray luminosity of normal galaxies, normalized 
at z=0, from Ghosh & White (2001). (right) Observed evolution in mean B-band luminosity for the HDF-N galaxies, with the linear fit to 
log L B versus z plotted as a solid line. The error bars shown were derived by the dispersion of luminosities in each redshift bin. 
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Fig. 3. — 0.5-2.0 keV and 2-10 keV flux distributions predicted for HDF-N galaxies. The lines and errors are deflned as in Figure 1. 



dominated by galaxies in this redshift range). ies expected to be detected as a function of flux in ultra- 

We proceed now to the question of the numbers of galax- 
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Flc;. 4. — The logN-logS for galaxies derived from the HDF-N estimates for the "Gaussian" (diamonds) and "Peak-M" evolution models in 
the 0.5-2.0 kcv (left) and 2.0-10.0 kcV (right) bandpasses. The solid lines give the logN-logS double power-lawr fit from Tozzi et al. (2001) 
based on the Chandra observation of the HDF-S (note that this logN-logS includes all contributions to number counts). The square points 
mark the source density determined in the 500 ks Chandra HDF-N observation (Brandt et al. 2001). Fx = 10^* erg cm~^ corresponds 
to ~ 24. 



deep Chandra, XMM-Newton, or future X-ray surveys. To 
address this, we derived the "logN-logS" number count 
distribution (i.e., the number of galaxies exceeding flux S 
as a function of flux) implied by our analysis, shown in 
Figure 4, using the HDF-N area sited above. The galaxy 
number counts have been corrected for spectroscopic com- 
pleteness using Figure 1 from Cohen ct al. (2000). This 
figure also shows the fits to the 0.5-2.0 keV and 2-10 
kcV HDF-S logN-logS distrubitions (Tozzi et al. 2001) 
(note that these curves are extrapolations below fluxes of 
~ 10~^^ and 10~^^ ergs cm~^ s~^ for the soft and hard 
bands, respectively). The 2-10 keV HDF-S logN-logS in 
the S — 10^^'"' — 5 X 10"^'* erg cm~^ s"""^ range is in 
good agreement with the HDF-N logN-logS curve given 
in Garmire et al (2001), where fluctuation analysis im- 
plies that the logN-logS curve flattens somewhat below 
S = 1.0 X 10^^"'' erg cm^^ s^^. The total (i.e., including 
AGN) source densities determined by Brandt et al. (2001) 
are marked which also implies a flattening of the X-ray 
logN-logS. Note that the results of the stacking analysis in 
Brandt et al. (2001) cannot be plotted on this diagram 
since by design the analysis was done around optically- 
detected galaxies without significant X-ray counterparts 
and hence the area is not well defined. Depending on the 
amount that the AGN logN-logS flattens, the total X-ray 
logN-logS may become dominated by normal galaxies at 
fluxes below logi^x = —17 18 erg cm~^ s~-^, particu- 
larly in the soft band. 

4. DISCUSSION 

We have estimated the X-ray fluxes of galaxies in the 
Hubble Deep Field using the known X-ray/B-band lumi- 
nosity correlation and taking evolution in the X-ray binary 
population of the galaxies into account (using the "Gaus- 
sian" and "Pcak-M" models in GWOl). Here we arc only 
attempting to estimate the LMXRB contribution to the 
galaxy fluxes, which should dominate the X-ray fluxes un- 
less significant starburst and/or AGN activity is present. 
The mean 2-10 keV X-ray flux that the Gaussian model 
predicts is ^ 8.2 x 10^^^ erg cm~^ s^^ which can be com- 
pared with ~ 1.7 X 10~^^ erg cm~^ s~^ obtained from a 
stacking analysis of a 500 ks Chandra observation of the 



HDF-N (converted to the 2-10 keV bandpass from a signal 
of 2.3 X 10~^^ erg cm~^ s~^ in the 0.5-8.0 keV bandpass 
assuming the same spectrum as discussed above) deter- 
mined by Brandt et al. (2001). The "Pcak-M" model 
predicts a mean X-ray flux a factor of ~ 3.5 lower than 
the "Gaussian" model. The "Gaussian" model is there- 
fore predicting a signal which is a factor of ^ 2 less than 
the signal observed, although we note that Brandt et al. 
(2001) selected optically-bright galaxies (a sample of 11) 
for their stacking analysis, and our estimate of the mean 
galaxy X-ray luminosity (prior to the application of any 
evolution) is evidently somewhat lower than that assumed 
by Brandt et al. (2001) but is nevertheless consistent with 
other estimates. In both the hard and soft bandpass, we 
are predicting that galaxies make up ^ 4% of the X-ray 
background down to a flux of ~ 10~^^ ergs cm~^ s~^. 

The best test of these models would be the direct de- 
tection of high-z galaxies in the X-ray bandpass, obviously 
with suflacient spectroscopic identification to rule out an 
AGN contribution to the X-ray flux. This would require 
an X-ray detection sensitivity on the order of 10~^^ erg 
cm^^ s^^, which would imply an exposure on the order 
of 10 Ms for Chandra (i.e., the effective exposure in the 
stacking analysis was 5 Ms), while even such a large ex- 
posure would be insufficient for XMM-Newton to reach 
these detection scnsitivites (due to XMM-Newton becom- 
ing background- limited at very large exposures). Opti- 
cal (and, to a lesser extent, IR) spectroscopic identifica- 
tions could possibly miss highly-obscured AGN or low- 
luminosity AGN which would produce enhanced X-ray 
fluxes relative to normal galaxies. A simple expectation is 
that the X-ray emission of galaxies would be extended (on 
spatial scales comparable to the optical emission), while 
the X-ray image of galaxies dominated by AGN would of 
course be unresolved. The half-light radius of galaxies in 
the HST Medium Deep Survey was typically ~ 0.1 - 1.0", 
with a median value around 0.6" (Ratnatunga, Grifiiths, 
& Ostrander 1999). Arcsecond resolution will therefore 
be necessary to detect any extension to galaxies. 

We thank the referee, Paulo Tozzi for useful comments 
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